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RATE OF REACTION BETWEEN MOLECULAR HYDROGEN 
AND MOLECULAR OXYGEN 
by Richard S. Brokaw 
Lewis Research Center 


SUMMARY 


The shock tube data of Jachimowski and Houghton have been rigorously analyzed to 

kn i 


obtain rate constants for the candidate initiation reactions H 2 + 0 2 H + HOg, 
k Q1 = 1.9xl0 13 exp (-24 100/T); H 2 + 0 2 -^-H 2 0 + O, k Q2 - 4. lxlO 13 exp(-25 400/T); 
and H 2 + 0 2 -22-OH + OH, k Q3 = 2.3xl0 13 exp (-25 200/T). (Rate constants in 

O 1 1 

cm mole sec , temperature in degrees kelvin. ) Reaction (01) is probably not the 
initiation process because the activation energy obtained is less than the endothermicity 
and because the derived rates greatly exceed values inferred in the literature from the 
reverse of reaction (01). Reactions (02) and (03) remain as possibilities, with reac- 
tion (02) slightly favored on the basis of steric and statistical considerations. The solu- 
tion of the differential equations is presented in detail to show how the kinetics of other 
ignition systems may be solved. 


INTRODUCTION 

When a mixture containing hydrogen and oxygen is subjected to a temperature and 
pressure pulse in a shock tube, small concentrations of atoms and free radicals are 
first formed from molecular hydrogen and oxygen. Reactions which have been suggested 
include (ref. 1) 

AH 298 /R - 28 700 K (01) 

AH 298 /R = 900 K (02) 

AH 298 /R = 9400 K 


H 2 + + H0 2 


H 2 + °2 


T)2 


h 2 o + O 


^03 

H 2 + P 2 -^0H,+ oh 


(03) 



The dissociation processes 


H 2 + M — 2H + M 


0 2 + M — 20 + M 


are sufficiently endothermic that they are unimportant except, possibly, at very high 
temperatures. 

After small concentrations of H, OH, or O are formed from reactions (01), (02), or 
(03), the atom and radical concentrations grow exponentially via the well-known branched- 
chain scheme 


OH + H 2 — L H 2 0 + H 


h + o 2 — 4-OH + O 


(I) 

(II) 


o + H 2 — iOH + H 


(m) 


At low temperatures and high pressures the chainbreaking reaction 

k 4 

H + 0 2 + M-4-H0 2 + M 


(IV) 


must also be considered. 

Recently Jachimowski and Houghton (ref. 2) studied the hydrogen- oxygen initiation 
process behind incident shocks. They analyzed their data to obtain rate constants as- 
suming initiation by reaction (01) or (03). Their analysis used an approximate and intui- 
tive formulation proposed to explain hydrogen- oxygen ignition delays (ref. 3). 

In this report the data of Jachimowski and Houghton are reexamined, using a rigor- 
ous formulation of the initiation and chain-branching kinetics. Rate constants are ob- 
tained by assuming initiation by reaction (01), (02), or (03). The values of rate con- 
stants indicate the more likely initiation processes. The solution of the differential 
equations is presented in detail to provide a guide as to how the kinetics of other similar 
ignition systems may be solved. 


THEORETICAL CONSIDERATIONS 

The differential equations governing the growth of radical concentrations during the 
induction period are as follows (ref. 3): 
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(1) 


= - (v 2 + V 4 )[H] + ^[OH] + ^ 3 [0] + ij 


d[OH] 

— = i/ 2 [H] - i/jtOH] + iUO] + 2i 3 

dt ^ 


d[0] 

dt 


= 


^[ 0 ]+i 2 


( 2 ) 

(3) 


Here [H], [OH], and [0]are the concentrations of hydrogen atoms, hydroxyl radicals, 
and oxygen atoms . Also, ^-^[Hg], “ ^2^2 ]’ ^3 = k 3^2^ V A. = and 

*! sk oi[ H 2][°2^ h = k 02f H 2 Jt°2 ^ i 3 s k 03t H 2^°2^ where t H 2^ and f°2J are the con " 

centrations of molecular hydrogen and oxygen, [M] is the total gas concentration, and 
the k's are the specific reaction rate constants for reactions (01) to (03) and (I) to (IV). 

During the induction period the concentrations of H, OH, and O build up rapidly, 
while the concentrations of H 2 and 0 2 are scarcely depleted. Hence, the v's and i’s 
in equations (1) to (3) may be taken as constants. 

The initiation rates ij, i 2 , and i 3 can be eliminated from the differential equations 
by introducing new variables C H s [H] + a H , Cq H s [OH] + aQ H , and Cq = [O] + a Q , 
where a.Q, aQ H , and aQ are constants. If these new variables are substituted into equa- 
tions (1) to (3), the initiation rates are eliminated by equating the sums of the constant 
terms to zero: 


("2 + "^H ' "l a OH " "3 a O + h = 0 

- i/ 2 a H + ^aQjj - i/ 3 a Q + 2i g = 0 

- ^ 2 a H + ^ 3 a 0 + i 2 = 0 
The new differential equations are 

dC H 

— — = ~(V 2 + )C H + ^Cqjj + I'gCQ 
d C/-\TT 

— — = ^ 2 C H - J'jCqjj + i/ 3 Cq 

dc o 

— - = U 2 C H. “ "3 C 0 


(4) 

(5) 

( 6 ) 

(V 

( 8 ) 

(9) 
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A particular solution to equations (7) to (9) is (refs. 3 and 4) 


C i = exp(Xt), i - H, OH, O (10) 

Substitution of equations (10) into equations (7) to (9) yields the relations 

~(i>2 + v 4 + A)A h + ^Aq h + ^Aq = 0 ( 11 ) 

^2 A H “ (^i + X ) A OH + ^3 A 0 = 0 ( 12 ) 

v 2 a h - (i^ + A)Aq = 0 (13) 


A nontrivial solution requires that the determinant of the coefficients of the A^ be zero, 
which leads to the cubic equation 
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A + (v-y + v 2 + ^3 + v 4 ) x + ^1^3 + v l v 4 + ~ v \ v 2^ v 2 " ^ 4 ) = 0 (14) 

In the chain-branching region (2^ > Va) Descartes’ rule of signs indicates that, if the 
roots of equation (14) are real, there are one positive and two negative roots. It can be 
shown that the roots are always real provided kj > 2kg, which is true for temperatures 
below about 2280 K, and hence true for the data of Jachimowski and Houghton (ref. 2), 
which are in the range 1100 to 1900 K. 

The general solution to equations (7) to (9) is 

Cj = A i j exp(Ajt) + Ag i exp(Agt) + Ag j exp(Agt) i = H, OH, O (15) 

where Aj, Ag, Ag are the roots of equation (14). We are specifically interested in the 
growth of hydroxyl concentration, which was the quantity monitored in reference 2. If 
we designate Aj as the positive root and Ag, Ag as negative roots, then 

[OH] = Aj q R exp(Ajt) + Ag QH exp(Agt) + Ag QH exp(Agt) - a QH 

~ A l,OH exp ^ A l^ ( 16 ) 

The approximation shown is a good one, except very early in the reaction, because a^ 
is small and the negative exponentials soon die out. 

We will now evaluate the coefficient A^ qjj. Initially, at t = 0, the hydroxyl radical 
concentration is zero, and equation (16) gives 

[° H ]t=0 = 0 = A l,OH + A 2, OH + A 3, OH " a OH 
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One can write similar expressions for hydrogen and oxygen atoms and then use any two 
of equations (11) to (13) to eliminate the and Aq in favor of Aq^. The same result 
can be obtained more easily by successive differentiation of equations (2) and (16): 

I = X l A l,OH + X 2 A 2,OH + X 3 A 3,OH = 2i 3 ( 18 

t=0 


d[OH] 

dt 


and 


d [OH] 2 2 2 

= X 1 A 1 , OH + X 2 A 2 , OH + X 3 A 3 , OH 


dt 


t=0 






=o 3 v dt 4=0 


- i/gij - 2 ^i 3 + v 3 1 2 

Equations (17) to (19) are solved to obtain 

X 2 X 3 a OH " 2 V X 2 + V + U 2 { 1 + V2 _ 2u l 1 3 


1, OH 


(^2 " x i )( x 3 _ 


(19) 


( 20 ) 


where the constant 


OH 


[2i/ 2 (ii + ^) + ^4(^2 + 2i 3)] 
vA2v 0 - v.) 


( 21 ) 


is obtained by solving equations (4) to (6). Expressions for Ag q H and A^ q h can be 


obtained by permitting the indices on Aj, Ag, and 


Ag in equation (20). 


Thus, equations (16), (20), and (21) describe the growth of hydroxyl concentration 
until such time as the effects of depletion of molecular hydrogen and oxygen are impor- 
tant or until the temperature rises due to atom and radical recombination processes. 


ANALYSIS OF EXPERIMENTAL INDUCTION TIMES 

The induction times reported by Jachimowski and Houghton (ref. 2) correspond to 
the time at which the hydroxyl concentration has risen to 10“ mole per cubic centi- 
meter. They also report experimental values of the growth constant. Thus, one can 
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obtain experimental values of qjj from equation (16) 

A 1 OH = 10_1 ° exp(-Ar) 


( 22 ) 


where A is the experimental growth constant and t is the induction time. Rate con- 
stants for initiation by reactions (01), (02), and (03) were obtained from equation (20), 
assuming that only one of the initiation reactions was occurring: 


101 = (X2 ' X ‘ )<X3 " Xl) 


2 i^2 A 2 A 3 


1-1 


vA2v 0 - v.) 


+ Vc\ 


(23) 


i02 = ii^] (i2 ' Xl,<X3 " Xl) 


(2^2 + 


1-1 


z/j(2i/ 2 - ^) 


+ Vc, 


(24) 


‘° 3 = [ 5 ^ (X2 ' Xl)(>3 ' Xl> 


"4 X 2 X 3 


1-1 


^( 2^2 ~ 


~ x 2 ~ X 3 ~ v \ 


(25) 


In these calculations the rate constants for reactions I, n, and III were taken from 
reference 5. The rate constants for reaction IV was taken from reference 6. And A^, 
A 2 , and A^ were obtained by solution of equation (14). 

These rate constants were least-squares fitted to the Arrhenius equation 

k = A exp/-— \ (26) 

\RT/ 

Results are summarized in table I, where the Arrhenius equation parameters and their 
standard deviations are presented. 


DISCUSSION OF RESULTS 

In this section the three candidate initiation reactions will be discussed in turn, with 
an indication as to which are the most likely initiation processes. 

The reaction 


H 2 + 0 2 - H + H0 2 


( 01 ) 


is a simple abstraction or two center reaction involving the breaking and formation of 
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one bond. Hence it is the most likely candidate based on steric considerations. How- 
ever, as Jachimowski and Houghton (ref. 2) have already observed, the activation en- 
ergy, 24 100 K, is substantially less than the endothermicity of 28 700 K. This amounts 
to more than three standard deviations in the activation energy and would seem to elimi- 
nate reaction (01), barring a temperature dependent systematic error in the experimen- 
tal data. 

Further, in reference 7 the rate of reaction (01) is reported to be 

k Q1 = 5.5X10 13 exp(-29 100/T) (27) 

in the range 290 to 800 K, with an uncertainty of a factor of 2.5. Equation (27) predicts 
a rate at 1100 K which is less than one-thirtieth of the rate at that temperature given by 
the constants in table I. Thus it is very unlikely that reaction (01) is the initiation pro- 
cess in these mixtures. 

The reaction 


H 2 + 0 2 - H 2 0 + O (02) 

might be termed a three center reaction inasmuch as new bonds are formed among three 
atoms. This would seem sterically more probable than the four center reaction 

H 2 + 0 2 — OH + OH (03) 

In addition, reaction (02) is less endothermic than reaction (03), although in either case 
the endothermicities are substantially smaller than the activation energies. Finally, 
the standard deviations in In k, In A, and E/R are slightly smaller for reaction (02) 
than for either reaction (01) or reaction (03). Initiation in these mixtures may be due to 
either reaction (02) or reaction (03) (or both), with reaction (02) the slightly more prob- 
able candidate. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 10, 1972, 

502-04. 
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TABLE I. - RATE CONSTANTS OF POSSIBLE HYDROGEN-OXYGEN 


INITIATION REACTIONS 




NASA-Langley, 1973 33 E-7195 
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